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Abstract 

We investigate the statistical properties of primordial black hole (PBH) 
formation in the very early Universe. We show that the high level of inho- 
mogeneity of the early Universe leads to the formation of the first generation 
PBHs. This causes later the appearance of a dust-like phase of the cosmologi- 
cal expansion. We discuss here a new mechanism for the second generation of 
PBH formation during the dust-like phase. This mechanism is based on the 
coagulation process. We demonstrate that the blue power spectrum of initial 
adiabatic perturbations after inflation leads to overproduction of primordial 
black holes with 10 9 g< M < 10 15 g if the power index is n > 1.2. 

1 Introduction 

Blue primordial power spectra of initial scalar perturbations have spectral index 
n > 1 and arise naturally in the modern inflationary scenario [1]. The observational 
upper limit on n is derived by normalizing the amplitude and spectral index n to 
the quadrupole scale. The COBE data [2] are compatible with a power spectrum of 
adiabatic perturbations P(k) cx k n with n = 1.2 ± 0.3at 68% confidential level (CL). 
At 95% CL the 4p observational limit on n corresponds to n < 1.65 [6]. It means 
that a direct extrapolation of the COBE data to extremely small scales even with the 
maximal possible value n = 1.5 (one-sigma upper limit) can give S rms great enough 
for the formation of a large (remarkable) number of mini black holes [3,4]. Limits 
on their production [5] can be used to re-evaluate the restrictions on the density 
perturbations spectrum including the spectral index n. This problem was discussed 
by Carr, Gilbert and Lidsey [3], and recently by Green and Liddle [4]. They found a 
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significantly strong restriction on the power spectrum index n < 1.3, rather smaller 
than n = 1.5 (one-sigma upper limit) from COBE-data. 

In our paper we demonstrate that due to the so called mechanism of induced 
formation a blue spectrum of the initial perturbations with 1.2 < n < 1.3 leads to 
the formation of a great number of PBHs with 10 9 g< M < 10 15 g 

Firstly, following [3] we assume that the PBHs of extremely low mass M in <C 10 9 g 
are formed directly in the early Universe from initial density fluctuations with power 
spectrum n > 1.2 and COBE-normalized amplitude at the moment to ~ ^-P'(j^")> 
where tpi and Mpi are the Planck parameters. We note that for blue spectra with 
n > 1.2 the variance of fluctuations af ms is great enough for the formation of a large 
number of mini black holes, which can produce a dust-like phase of the expansion of 
the Universe at times t 3> t . 

Secondly, we are interested in the situation where PBHs with 10 9 < M < 10 15 g 
are formed from PBHs with M in <C 10 9 g due to hierarchical clustering of fluctuations 
in the dust-like background of the primordial low mass black holes (M = M in ). We 
demonstrate that during a PBH dust era the fraction of the matter in the Universe 
accumulated into PBHs with 10 9 < M < 10 15 g is great enough to produce the 
remarkable abundance of the massive PBH in the indicated interval of masses. These 
black holes will evaporate up to the present time and will be able to change the 
observational restrictions [5]. 

Finally, we shall demonstrate that for a blue power spectrum of initial adiabatic 
perturbations there is a limit on the power index: n < 1.23. 

The plan of this paper is as follows. In Section 2 we briefly discuss the conditions 
for PBH formation in the early Universe with the equation of state p = 7P and 
7 = 1/3. In Section 3 we investigate the possibility to have a dust-like phase of the 
expansion of the Universe due to PBH creation with extremely low masses M in . In 
Section 4 we discuss the conditions for formation of PBHs by direct collapse of the 
initial perturbations. In Section 5 we investigate a more efficient mechanism of PBH 
formation during the dust-like phase of the expansion of the Universe which is based 
on the coagulation approach. In Section 6 we summarize the main results of the 
investigation. 

2 Formation of the PBHs with the extremely low 
mass. 

In this section we investigate primordial low-mass black holes formation from density 
perturbations directly at the epoch of the expansion of the Universe with the hard 
equation of state p = jp. We assume that just after the period of reheating, which 
is a natural part of the modern theory of inflation, the equation of state corresponds 
to 7 = 1/3 and the spectrum of the initial perturbations is blue: P(k) ~ k n ; n < 1.5 
according to COBE data (one-sigma upper limit) [1]. Following [6] we assume that 
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initial density perturbations produced during the inflation period are Gaussian and 
the probability distribution of the smoothed density field P(5(M)) is given by 



P(5(M))dS(M) = -= — - exp 

V2ira rms (M) 



5 2 {M) 



d5{M). (1) 
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Here a(M) is the mass variance evaluated at the horizon crossing; 

1 POO 

a 2 (M) = / dkk 2 P(k)e~ k2R2 ^ (2) 

Z7T Jo 

where P(k) = A 2 k n is the power spectrum with amplitude A 2 and spectral index n; 
R(M) is the filtering scale which is linked with the mass of fluctuation Mh = ^-pH~ 3 
at the horizon scale R ~ M 1//3 for 7 = of the equation of state and R ~ M 1//2 for 
the radiation-dominated Universe. In the 7 = 1/3 Universe the mass fraction of 
black holes /5(M), which are formed directly by the collapse of initial perturbations 
at the horizon scale can be obtain from formula (1) using Carr-Hawking criterium 
I < S(M ~ M H ) < 1 (see [7]): 



P(M)= f 1 d5P[5(M)]~a rms (M)exp 

Jl/3 



l8a 2 rms (M) 



(3) 



Note, that there are a number of observational limits on a PBH fraction in the 
very early Universe, covering various mass ranges of its mass spectrum (3). For the 
future investigation we use / 9(M)-limits from nucleosynthesis and PBH evaporation 
at present time [5]. As one can see from (3), mass fraction of PBH f3(M) = p p bh ^ Min ^ 
is very sensitive to the amplitude and spectral index n of the initial density pertur- 
bations. Here p p bh,in{Mi n ) is the mass density of black holes and p 7 is the density of 
the ultrarelativistic background at the moment of PBH creation. 

For example, if a rms (M) oc 1CT 5 and n — 1 for Harrison=Zeldovich spectrum 
normalized to COBE-data, then a probability of PBH formation is extremely low: 
(3 ~ 10~ 6 ' 10 ~ 0, and primordial black holes practically cannot exist in the whole 
mass interval. However, the situation changes drastically for blue power spectra 
n > 1. It is worth to note , that for 7 = 1/3 Universe the variance of fluctuations at 
M scale has a form [4] : 

/ M \ 2 

a 2 rms (M) = a 2 ms (M eq ) \=-J ; M < M eq , (4) 

where M eq is the standard horizon mass at the moment of matter-radiation equality. 
For M > M eq up to COBE mass scales, the variance is 

n-1 

a 2 rms (M) = a 2 rms (M eq ) [=-J ; M > M eq . (5) 
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For the future investigation we use Eq.(4) and Eq.(5) and the fit of the four-year 
COBE data for the amplitude of spectra at the scale K = a H from [8]: 



A = V2 n{avH Q )-^5 H , 

5 H = 1.9 x 1(T 5 exp[1.01(l - n)}, 

where H and a are the Hubble constant and scale factor at the present time. The 
distribution of function a rms (M) for different M, is plotted on the power index -mass 
diagram (Fig.l). Further we use the modified power index a = v ^ and normalization 
a rms (M) from Eq.(4)-Eq.(5) on the COBE mass scale M C obe, which corresponds to 
K = a H : 

Let us introduce the new function a (M,a rms ): 



lg 


& rms 




J?rms{McOBE) _ 
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a{M,a rms )= — -. r-, (7) 



where a rms = const - is the parameter. 

Note, that for any kind of blue power spectrum the power index a = a* = const. Let 
us consider the equation a (M, a rms ) = a*. For different values of the parameter a rrns 
this equation gives us the function M = M*(a rms ), which corresponds to the a rms (M) 
from Eq.(6) at a = a». Using this method, we would like to focus attention on the 
thin dashed, dotted and solid lines of Fig.l, which are concentrated at the range of 
a < 0.09 (see Fig.l). These lines correspond to the series of function a(M, a rms ) 
from Eq.(7) with the different values of a rms . The lowest dashed line corresponds to 
&rms = 0.020, the next line (dotted) corresponds to a rms = 0.025. 
Three next solid lines represent a rms = 0.030; 0.035 and 0.040 correspondingly. Two 
next lines (short dashed and long dashed) correspond to a rms = 0.045 and a rms = 
0.050. As in can be seen from Eq. (3) and Eq.(7) for this kind of function a(M, a rms ) 
the value of parameter a rms can be normalized on the fraction (3 of PBH (see Eq.(3)). 
Using this normalization, we plotted in Fig.l the thick solid line at the interval 14 < 
lg M/Mpi < 20 which corresponds to the well known limits of PBH's probabilities 
[5]. ' 

There is no observational restriction at M < 10 9 g, and we can use only one limit 
Qpbh < 1, where fl p bh is the ratio of the black hole and critical densities at present 
time, for M > 10 15 g . 

Note, that the increasing of the variance of fluctuations at M — > Mpi is the 
natural behavior of the blue power spectrum. There is not such a behavior for 
Harrison-Zeldovich and tilted (n < 1) power spectra. The growth of a rms at small 
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scales of perturbations is a cause of the producing low-mass PBHs (M <C 10 9 g but 
M ^> Mpi) which evaporate before the period of the cosmo logical nucleosynthesis. 
However, these low-mass PBHs can produce the dust-like phase of the expansion of 
the Universe and change the condition of formation the next generation of PBHs 
with 10 9 g< M < 10 15 g. Thus, the high density of the low-mass PBHs (we use an 
index "m" for this first generation of PBHs) can induce the overproduction of the 
PBHs with 10 9 g< M < 10 15 g (the second generation) at the dust-like phase of the 
expansion. 

In the next section we will discuss this possibility in some detail. 



3 Dust-like phase of formation of the PBHs. 

Let us describe the condition which produces the PBH dust-like phase of the expan- 
sion of the Universe. We denote by t in the time of formation of black holes of the 
first generation in the radiation dominated Universe. 

Let (3(M in ) be a fraction of matter to the Universe collapsed into the first gen- 
eration PBHs with mass M p bh = It is well known [4, 6] that the period of 
formation, tj n , and the mass of PBH M = Mi n are linked by the equation: 

M m = 1 ^M H = 1 -^( t -^)M Ph (8) 



g 1 / 2 \tpiJ 

where gf is the effective number of massless degrees of freedom at t — t in . After 
primordial black holes are formed, the mass density of its first generation p p bh,in 
begins to decrease as p p bh,in oc a -3 , where a is a scale factor. The mass density of 
ultrarelativistic matter decreases more rapidly: p 7 oc a~ 4 . In this case at 

t > t eq , t eq ~ t in {M in )(3- 2 {M in ) (9) 

the cosmological expansion is determined by PBHs with M = Mj„. If the time of 
evaporation of these black holes 

1.2-10 4 /M in \ 3 , , 

*■«, = —— hrr 1 -tpi, (10) 



K \Mpi 

is much greater than the time t eq , then the duration of PBH dust-like phase is At = 
r ev (M in ) ^> t eq (M in ), where n is the parameter depending on the number of particles 
species which can be emitted (see [4])). 

It is worth to note that at the moment t = ti n the space distribution of PBHs 
with M = M in is modulated by the long wave perturbations in the ultrarelativis- 
tic background. This spectrum of adiabatic fluctuations is preserved in the PBH's 
distribution up to the dust-like phase. 
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At t > t eq density fluctuations in the non-relativistic "gas" of PBH growth with 
time and the variance of perturbations in the PBH's distribution at t >> t eq is: 

( t \ 2/3 

a rms {M, t) = a rms {M eq )T{M, M eq ) I— J (11) 
where M is the mass of perturbations, and 




2a 
3 



at M < M eq 



T(M,M eq ) = { ^ m+a) ~ (12) 



3 at M > M 



eq , 



eq 



is the transfer function of perturbations, and M eq is the mass of PBHs in the volume 
V ~ t 3 : 

M eq ~ (3-\M in )M m (13) 

As it is seen from (11) - (13), the spectrum a rms (M,t) is practically flat (a <C 1) 
when M < M eq and there is a turn over point at M > M eq . This spectrum arises only 
under conditions M in <C M eq and /3(M in ) <C 1, when the hydrodynamic approach is 
correct. Actually (see Fig.l) for the blue power spectrum with n = 1.23 the level of 
primordial fluctuations at small mass scales is comparatively high. For example, if 
M in w lg, then a rms (M in = lg) w 0.1; /3(M in = lg) w 4.4- 10~ 4 and M e9 w 5-10 7 M m . 
In this case the spectrum of perturbations in the black holes background is similar 
to well known spectral distribution of fluctuations in the standard Cold Dark Matter 
(CDM) Universe model. 

This analogy has a deep reason. In the standard CDM theory of the formation of 
the large scale structure of the Universe the typical scale of the mass density pertur- 
bations has a turn over point at the scale of the equality, which is close to the typical 
scale of clusters of galaxies (M**). At the scale M < Mf the dynamics of cluster- 
isation of perturbations is very complicated because of nonlinear transformation of 
spectra which was investigated by Doroshkevich et al. [9]. 

According to [9], the hierarchical clusterisation of non-linear configuration takes 
place at all scales of perturbations. During this process the typical configurations 
of low mass (mini pancakes) are clustering step by step into massive pancakes up to 
the typical scale M Mf q . In the standard CDM Universe clusters of galaxies are 
formed during a very short time corresponding to a redshift z, < z < 1. It means 
that nonlinear transformation of the spectrum of perturbations has a cutoff caused 
by the limitation of the age of the Universe. 

For the PBH's dust-like phase at r ev ^> t in the situation changes drastically. The 
duration of this phase covers not only the period of the formation of structures with 
M ~ M eq , but also the period of the formation of structures with M = M max ^> M eq . 
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In this case objects with M pa M eq at the moment of evaporation r ev are virialized 
and form the high massive objects with M = M max (see (11) - (12)): 

i 

M max ~ M eq a^ (M m ) {^^j ^ ■ (14) 

In addition, there are two very important differences between the rate of cluster- 
isation during PBH dust-like phase and the standard CDM model of the structure 
formation. Firstly, if (5{M in ) is not especially small, then the scale of maximum of 
perturbations M eq contains N < 10 2 mini PBH with M = Mj n , and for the descrip- 
tion of the rate of clusterisation one needs to use the kinetic approach. This situation 
takes place, for example, if initial power index is n = 1.26 and M in =lg (see Fig.l). 
As it is seen from Fig.l, a rms (M in = lg) = 0.3; (5{M in ) ~ 0.16 and N ~ 40. 

Secondly, in the context of the CDM scenario the background particles are colli- 
sionless and "interact" each other and with the baryonic matter only through gravita- 
tion. However, PBHs with M = Mi n can collide especially inside the quasi stationary 
stable gravitational configurations (PBH's clusters - "gas"), which are formed during 
the dust-like phase. 

These two last reasons are linked with a very important question - can the first 
generation primordial black holes produce the next generation of the PBH during the 
dust-like phase or they have to disappear because of Hawking's evaporation without 
any observational consequences? 

The answer this question depends very essentially on the probability of PBH 
formation at the dust-like phase, which we describe in the next section. 



4 The probability of the formation of the second 
generation of PBHs. The spherical symmetry 
approximation. 

The problem of PBH formation during the dust-like phase was pointed out in the 
pioneering articles by Polnarev and Khlopov [10,11]. According to [10] the fraction of 
matter going into PBH is determined by the probability that the overdense region is 
sufficiently spherical symmetric and homogeneous. This region collapses into a black 
hole directly. The fraction of PBH which is formed due to this mechanism is given 
by [10]: 

P(M bh ) ~ 2 • 10- 2 [a rms (M,t,)] 13 / 2 , (15) 

where a rms (M, to) is a variance of the density perturbations at the moment td when 
the dust-like phase starts. 

Below we briefly discuss the initial conditions leading to Eq.(15). As we mentioned 
above, if the hydrodynamic approach is applicable, we can describe the statistic 
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properties of the Gaussian field of the density perturbations according to Bardeen et 
al. [12]. First of all, we can express the Gaussian point process entirely in terms of 
the density perturbations field 5(f) and its derivatives. In the vicinity of the point of 
a maximum of 5 field f=fl we can use the Taylor series expansion: 

*(r) = S(f*) + W KM ~ r*l (rj - r*) (16) 

where Ay is the second derivative tensor of the field 5 at f = r* (so called defor- 
mation tensor). We transform the coordinate system along the principle axes of the 
deformation tensor and determine three its eigenvalues Ai, A2 and A3. Doroshkevich 
[13] has shown that the joint probability density for Ai, A2 and A3 has a form: 

P(Ai,A 2 ,A 3 )= (17) 

3/i 2 - 7.5/i 2 \ 



£>(Ai- A 2 )(Ai- A 3 )(A 2 - A 3 )exp 



<7in 



where D = K^ a z 7 , &in is the variance of 5(f) at the moment t^, /xi = Yh=i \ an d 

in 

/i 2 = A1A2 + A 2 A 3 + A1A3. Following [12] we define a three-axial ellipsoidal surface 
with semiaxes 

-2(5(0) -Sir)) 1 ' 2 



A ? ; 



(18) 



and characterize the asymmetry of 5(r) field at the vicinity of the point of maxima 
(r* = 0) by the parameters: 

As the result, the distribution of 5(r) field around a point of a maximum in the 
spheric system of coordinates is [12] 

2 

5(r) ~ 5(0) - xa 2 r -[l + D(e,p)]; r -> 0, (20) 

where 

D(e,p) = 3e[l - sin 2 ^(l + sinV)] +p[l ~ 3 sin 2 ^ cos V], 

x = -^l; al = ^jdkk^P(k) (21) 

P(k) is the power spectrum of the density perturbations. 

The main idea which was used by Polnarev and Khlopov [10,11] for the determina- 
tion of the probability of PBH formation during the dust-like phase of the expansion is 
based on the hypothesis that three eigenvalue of the deformation tensor are extremely 
closed to each other: Ai — A 2 — Ai — A3 ~ a in . Thus the parameters e and p are going 
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to zero, and the distribution of the density field around the point of the maximum 
is spherical symmetric with a very good accuracy. One can obtain the probability 
of such realization of the random Gaussian process by integration of Eq.(17) over 
A2 and A3. After integration we have /3(Mbh) oc <rf n [10,11]. An additional factor 
in Eq.(15) is connected with the condition of the homogeneous collapse of the 
configuration according to the Tolman's solution [10,11]. Taking into account these 
conditions we would like to introduce two additional criteria which are important for 
the determination of the probability of PBH formation during the dust-like phase. 

First, the collapse of a configuration into a black hole is definitely a non-local 
process and a criterium of the field distribution 5 around a point of a maximum 
should be non-local also. Usually, this non-locality is possible to imitate using the 
filtering procedure: 



where Ri are different scales of filtering. Let us focus our attention on the spectrum 
of perturbations in the PBHs (or another "dust" particles, as was investigated in 
[10,11]), using Eq.(12). There exists a low-massive plateau at M < M eq , where 
a rms (M < M eq ) is practically constant for a <C 1. 

In this interval we use the small scale of filtering R s . At the same time, at the 
scale M > M eq in terms of Eq.(12) we will use the large scale filtering Ri. Following 
[10,11], we assume that after filtering (R = R t ) the condition of the spheric symmetry 
is valid in the vicinity of some point of a maximum of field S(f, Ri). The criterium of 
the homogeneity of distribution 5(r, Ri) is definitely correct for the .^-filtering only. 

At the small scale R ~ R s the situation is more complicated. Following [12], 
it is easy to show that the concentration of the peaks which is obtained by the 
filtering R ~ R s is much higher than the the concentration of peaks (not especially 
high) in 5{r, i?;)-field. The correlation function of the "smooth scale filtered" peaks 
has a typical correlation radius r ~ R s . If r ^> R s , then the character of the space 
distribution of the small scale peaks is Poissonian. It means that inside a large filtered 
peak which is spherical symmetric and homogeneous in average, there are N ^> 1 
small peaks which have a Maxwellian velocity distribution [12]. This configuration 
(a large scale smoothed peak and a lot of small scale peaks inside it) has properties 
to form not a black hole but a virialized cluster. However, the formation of a cluster 
is the natural part of the hierarchical clustering scenario without any specific and 
low-probability conditions of symmetry and homogeneity of the initial configuration. 

In the CDM model these clusters are formed as a result of the interaction of 
pancakes (see [9]). It is the first reason why the probability f3(M) from Eq.(15) 
overestimates the number of PBHs which can be formed during the dust-like phase 
of the expansion of the Universe. 

The second reason is the following. The field of peaks which is smoothed by Ri 
filter has a peak-peak correlation on a typical scale r > R[. It means that a local 




(22) 
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criterium of spheric symmetry, which does not take into account the influence of 
neighboring peaks, has a limited range of application in the especially small vicinity 
of a point of a maximum. We demonstrate this situation in Fig. 2. As it is seen 
from this figure, the influence of the neighboring peaks can practically destroy the 
spherical symmetry of an initial peak. This problem of an "outside interaction" is 
the second reason to decreasing of / 9(M)-probability. 



5 Induced (coagulation) mechanism of the PBH 
formation 

In our model the dust-like phase is linked with the first generation of the PBHs, and 
we can introduce the following new mechanism of the PBH's formation, based on the 
typical realization of overdense configurations - like pancakes. We don't need a special 
spheric symmetry and homogeneity of the configuration because the formation of 
clusters is the result of interactions between pancakes similar to the cluster formation 
in the standard CDM model. 

Using Press-Schechter formalism we can estimate the moment of formation of the 
PBH clusters for the spectrum (11) 

w a H uw ) <<UM '" ) ' (23) 

where a cr = 1.69. The majority of these clusters are practically virialized at t ^> 
t(M eq ). Let us estimate the characteristic time scale when the PBHs inside a clus- 
ter can collide. The cross-section of the interaction of PBHs is a pbh = 7rr|, where 
r g = 2GM in is the Schwarzschild radius, and the variance of the velocity is v% ir ~ 
GM eq / R eq , where R eq is the radius of a cluster. In this case the time scale of the 
"hard core" collision is: 

-j {M eq \ 2 ( a cr \ 7/2 

r co i ^ {(T P bhn pb hV V ir) = t ( — J -—— . (24) 

\M in J \a rms (M in ) J 

Note, that the r coi is much greater than the time of the velocity relaxation t v inside 
a cluster: 

r, = r co ^ r ln- 1 (^)«r coi . (25) 

It means that at t > t co \ the interaction between black holes of mass M = M in 
inside a cluster can be described using the coagulation approximation [14-16]. Let 
n(m, t)dm be the number of the density of black holes in the mass range m -r- m + dm 
at the moment t. For mass distribution function n(m, t) we can write the coagulation 
equation [14,15]: 

dn 1 f m 

— — - / a c (m' ,m — m! ,t)n(m' ,t)n(m — m! ,t)dm— (26) 
ut 2 jo 
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-n(m,t) / a c (m,m' ,t)n(m' ,t)dm'; 
Jo 



where a c = {a p bh{m, m')v), and ( ) denotes the averaging over the peculiar velocity v. 
After averaging over the peculiar velocity we have a = a p i,h(m, m')v V i r . From Eq.(25) 
we can obtain the following expression: 

dM tot r°° , ,dn 



dt 



[°° , ,0n 
I dm ■ m —— = (J: 
Jo dt 



dN f°° ,dn 
~dt 



= /„ *<? = < 27 > 

1 roo roo 

— — - / dm' / dma c (m,m')n(m.t)n(m ,t), 

2 Jo Jo 

where M tot is the total mass inside a cluster and iV is the total number of the struc- 
tured elements in the cluster. We would like to emphasize that the rate of reaction 
a c does not depend on the time t. This approximation is practically the same that 
was discussed by Trubnikov [14], and Silk and White [15]. In this case we can use 
the scaling invariance of the function a c (m,m'): 

a c (km, km') = k h 'a c (m,m')] v — 2, (28) 

and describe the distribution function n(m, t) as: 

^^"■•^(-^l (29) 

where m*(i) is the typical scale of mass during the process of coagulation of the 
PBHs inside the cluster, and F is some function which depends on (~ m "| t ) ) ■ This 
new function has the following property: F — > for m/m* — > oo. From 
Eq.(25)-(27) for v = 2 one can find the equations for m*(i) and N(t): 

d 4- = -C; ^ = Dml (30) 
dt ' dt y J 

where C and D are the constants of normalization. 
As it is seen from (29) N decreases with time as: 

N = Nq — Ct, (31) 

where N is the initial number of points at t — 0. In our case N Q = M eq /M in . For 
the determination of constant C we can use the initial condition of the coagulation 
process: 

n(m, t = 0) = N 5{m - M in ), (32) 
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where 5 is the Dirac ^-function. Substituting it in Eq.(32) in Eq.(26) we have 

^ = -\a c (M m )N 2 = -C. (33) 

Let us come back to Eq.(29). As one can find from this equation, mass of black holes 
of the second generation which are formed inside a cluster increases with time as 

M - = T^M-Df < 34 > 

and formally going to infinity at t = (M in D)~ l . Actually at t = jV| ^~ 1 the whole 
mass of a cluster is concentrated in one massive black hole. But in this limit the 
coagulation approximation with virialized velocity field is inapplicable. However, the 
rate of the coagulation can be used for qualitative estimation of the mass of a central 
black hole inside a cluster. At t ~ (iVo — 1)/C it is equal: 

M bh ^ ^Zdnq =V'M eq - ii< 1, (35) 
1 c 

and 

D ~ (2M m r col )- 1 (l - ^) . (36) 

We would like to emphasize that the rate of coagulation (31)-(32) is only the low 
limit of the clusterisation of the first generation of the PBHs into massive primordial 
black holes of the second generation. It is caused by the taking into account only 
"hard-core" collisions of black holes inside a cluster. In a general case the formation 
of pairs is important and the rate of coagulation increases^. Thus if r col <C r ev (M in ), 
then the PBHs form a massive black hole with M 3> M in inside a cluster. In Fig.l 
we plotted the thick solid line with diamonds which corresponds to the equality 
T C oi = T ev (Mi n ). If initial amplitudes of perturbations at the range Mi n < 10 14 g are 
a rms (M in ) > 0.1 then the first generation of the PBHs produces the dust-like phase 
and induce the formation of the second generation of black holes. Two thick lines with 
asterisks correspond to the condition M eq < 10 9 g (left line) and M eq < 10 15 g (right 
line). In the first case Mj„ has to be less than 10 5 Mp/. At the second case we have 
M in < 10 7 Mp h In both cases the critical level of the power index n is n = 1.23. What 
will happen if this condition is broken? Obviously the first generation of PBHs forms 
the second generation of PBHs with 10 9 g< M < 10 15 g and /5(10 9 g< M < 10 15 g) < 1. 
We can estimate this factor (3 using the Press-Schechter formula. For the clusters 
with M Pd M eq the concentration is the following 

,,,,,,„ 4 tw M V" + °' /3 



1 We describe this situation in a separate paper. 



12 



x exp 



M v 2(l+a)/3 - 



M max {t) i 



dM, (37) 



where M max is determined by Eq.(14), and p is the mean density of the PBHs during 
the dust-like phase. This density of the clusters with masses M eq < M < M max 
is Pd ~ ~Ppbh- K each cluster of mass M eq contains a massive black hole with the 
mass Mbh = pM eq , p < 1, then the probability of the second generation black hole 
formation is 

PHh(M bh ) ~ ii > (38) 
where (3^ is the observational limitQ. 

6 Conclusions 

We have presented here the new mechanism of the formation of the primordial black 
holes in the very early Universe after inflation. We have shown that a high level of 
adiabatic perturbations in the early Universe leads to producing the first generation of 
the PBHs. Further these PBHs form the specific phase of the cosmological expansion 
- so called dust-like phase. The hierarchical clustering of perturbations in a "gas" of 
the PBHs is the source of formation of the next generation of the black holes inside 
the clusters. We have presented the new mechanism of the second generation of the 
PBHs formation which is based on the coagulation mechanism. We have shown that 
the initial inhomogeneity of the Universe "switch on" the mechanism of the induced 
PBH formation and leads to their over-producing if the power index of perturbations 
is n > 1.2. We note that the formation of PBH from the "point-like" collisionless 
particles and from PBH of first generation has only one difference. Using Eq.(24) 
one can see, that t co i oc M~ 2 and goes to infinity, if M = m — > , where m is the 
mass of particles. In this case the coagulation mechanism of the PBH formation is 
inapplicable and the probability of the PBH formation is extremely low. 
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2 We would like to note that /i is remarkable less than unity, for example due to emission of the 
gravitational waves [17] 
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Figure 1: The diagram: the power spectrum index - mass. The tilted dashed, dotted, 
thin solid and so on lines correspond to the different values of the parameter a rms 
of the function a(M,a rms ), see formula (7). For seven of these lines in the bottom 
part of the diagram, from the bottom to the top the amplitude of perturbations a rms 
changes from 2 • 1CT 2 to 5 ■ 10~ 2 with the step A = 5 ■ 1CT 3 . The next series of the 
tilted lines (in the middle part of the diagram) corresponds to the levels of amplitude 
&rms — 0.1 -r 0.4 with A = 0.1. The last solid tilted line corresponds to the amplitude 
&rms = 1- Four horizontal lines from the top to the bottom correspond to n — 1.5 
(COBE-restrictions on the spectrum power index); n = 1.3; 1.26; 1.23. The solid line 
with asterisks inside the range 14 < lg M/M Pi < 20 corresponds to the observational 
restrictions [5]. The solid line with diamonds corresponds to the condition r col = r ev 
(see Section 5). Two solid lines with asterisks in the range 5 < lg M/Mpi < 14 
correspond to the condition Mbh < 10 9 g and Mbh < 10 15 g correspondingly (see 
text). The vertical line M = 10 14 Mpi is the lowest limit of mass for which there exist 
observational restrictions, based on the Hawking's process. See the text for other 
details. 
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Figure 2: The Distribution of the density field in the vicinity of the point of maxi- 
mum for the typical realization of the two-dimensional Gaussian field with Harrison- 
Zeldovich spectrum of perturbations. The origin of the system of coordinates corre- 
sponds to the point of maximum of 3ay ms -peak. The lines represent isodensity levels, 
x and y coordinates are normalized to the filtering radius Ri. There is a cluster of 
three peaks: 3.3a rms , 2.6a rms and 2.2a rms . Each member of the cluster determines 
the curvature and symmetry of the isodensity lines in their vicinities at the slice above 
« = 2-f-2.5 (see text). 
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